Impurity-induced bound states in iron-based superconductors with s-wave cos kx ■ cos k. 

pairing symmetry 
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Using both the self-consistent Bogoliubov-de Gennes formulation and non-self-consistent T-matrix 
approach, we perform a comprehensive investigation of the in-gap bound states induced by a localized 
single impurity in iron-based superconductors. We focus on studying signatures associated with the 
unconventional sign-changed s-wave pairing symmetry. For a non-magnetic impurity, we find that 
there are two in-gap bounds, symmetric with respect to zero energy, only in the sign-changed s-wave 
pairing state, not in the sign-unchanged s-wave state, due to the existence of non-trivial Andreev 
bound states caused by the sign change. For a magnetic impurity, we find that due to the breakdown 
of the local time-reversal symmetry, there exist only bound state solutions (with orbital degeneracy) 
carrying one of the electron-spin polarizations around the impurity. As increasing the scattering 
strength, the system undergoes a quantum phase transition (level crossing) from a spin-unpolarized 
ground state to a spin-polarized one. While the results for the magnetic impurity are qualitatively 
similar in both the sign-changed and sign-unchanged s-wave superconducting (SC) states, the bound 
states in the first case are more robust and there is no n phase shift of the SC gap near the impurity 
in the strong scattering regime. 

PACS numbers: 71.10.Fd, 74.20.-z, 74.25.Jb 



I. INTRODUCTION 

Iron-based superconductors, a newly discovered 
family of superconductors with layered conducting 
planes r''^'^'^ have stimulated enormously theoretical 
and experimental studies in the condensed matter com- 
munity. This is not only because of the Tc can be 
as high as 55K, but is also due to its striking simi- 
larity to the high-Tc cupratcs. One evident observa- 
tion is that the undoped iron-based superconductors, as 
cuprates, exhibit antiferromagnetic order (though still 
metallic) and the superconductivity occurs only when 
electrons or holes are sufficiently doped into the Fe lay- 
ers. Due to the close proximity to the magnetism, the 
mechanism of the superconductivity is expected to be 
unconventional^'^'^'iaaiiis^ii 

As an important hint to unveil the mechanism of the 
superconductivity, the determination of an explicit super- 
conducting (SC) gap structure is indeed essential. The 
basic electronic band structures of the iron-based su- 
perconductors have been predicted by the first-principle 
calculationsi^ii^ that there are two hole Fermi pockets 
around F point and two electron Fermi pockets around 
M point in the first Brillouin zone (FBZ) with all five 
orbitals of an iron atom getting involved. As a result of 
this complexity, undoubtedly, nailing down right SC gap 
structure is still one of the most challenging issues in this 
rapidly growing field. 

On the theoretical side, there have been many pro- 
posals for the possible pairing symmetries of iron pnic- 
tides, including nodeless or nodal SC order parametersJ^ 
Majority of studies so far, either from weak-coupling or 
strong-coupling approach, suggest an extended s-wave 
state (so called s± state), where the relative sign of SC 
order parameters changes between the hole and electron 



pockets (see Fig.[T|)iii^i2iii^ In particular, in a recent paper 
done by Sco et al.^ it gives, in the proposed two-orbital 
exchange coupling model ji^i an explicit cos fc^; - cos ky form 
of the pairing symmetry in momentum space, as long as 
two general conditions are satisfied: (i) the next nearest- 
neighbor (NNN) superexchange coupling J2 dominates, 
and (ii) Fermi pockets are small near the aforemen- 
tioned spots in the FBZ. Furthermore, this simple form, 
cos kx-cos ky, is in good agreement with the SC gaps mea- 
sured by the angle-resolved photo-emission spectroscopy 
(ARPES)i^i2ii^ Although this model does overestimate 
its insulating behavior near the undoped regime, this de- 
ficiency is irrelevant to determination for the properties 
of the SC state. We will focus on this pairing state (for 
simplicity, s± state) within this model throughout this 
paper. 

On the experimental side, however, there is still no cen- 
sus on the gap structure of the SC state. For instance, the 
ARPES results support a fully-gapped SC statcr^i^ii^ 
consistent with point-contact Andreev spectroscopy^^ 
and magnetic resonances measured by neutron scatter- 
ing experiments r^i^^i^ while some penetration depth 
experimentsr2i2^i^ NMR experiments^>^i^i2^ and 
other experiments^ seemingly contradict with former in- 
terpretation. Even though it has been argued that s± 
state could partially reconcile these difficulties f^^i^^i^i^^ 
it is still far from this SC order parameter being settled 
down. This, again, reflects an urgent call for a practical 
way to detect the SC gap structure, especially with sen- 
sitivity of measuring the sign change in the internal SC 
phase. 

Despite of several theoretical works proposing vari- 
ous ways to detect the phase structure of the SC order 
parameter achieving experimental realization 
reliably remains challenging. Since disorder is an intrin- 
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sic property in doped superconductors, a comprehensive 
study of impurity effects can also help indirectly probing 
SC order parameters. Here we propose tunneling mea- 
surements of impurity-induced states as a probe which 
is sensitive to the internal phase of the gap function be- 
tween electron and hole Fermi pockets. Such sort of mea- 
surements has been proved invaluable in determining the 
nodal d-wave pairing symmetry of the high- Tc cuprates4i 
In addition, the study of the impurity effects can also pro- 
vide useful information about the SC gap structure and 
even uncover competing orders j ^^i'^^'^'^i'^^'^^i^"'^"'^ 

In this paper, we perform a detailed investigation of the 
impurity-induced bound states in iron-based supercon- 
ductors within a two-orbital exchange coupling model. 
By using both the self-consistent Bogoliubov-de Gennes 
(BdG) formulation and non-self-consistent T-matrix ap- 
proach we find the following general features, (i) For the 
non- magnetic (intra-orbital) impurity potential, there ex- 
ist two in-gap bounds, symmetric with respect to zero en- 
ergy, only in the s-wave cos ■ cos ky pairing state, not 
in the sign-unchanged s-wave state. The basic physics 
of this result stems from the emergence of non-trivial 
Andreev bound state within the SC gap due to the im- 
purity scattering that destroys any unconventional (as 
opposed to the usual s-wave) phase assignment, (ii) For 
the magnetic (intra-orbital) impurity potential, due to 
the breakdown of the local time-reversal symmetry, there 
exist only bound state solutions (with orbital degener- 
acy) for one of the electron-spin polarizations around the 
impurity. As increasing the scattering strength, the sys- 
tem undergoes a quantum phase transition (level cross- 
ing) from spin-unpolarized ground state to spin-polarized 
one. Unlike the case with the non-magnetic impurity, the 
results in the s-wave cos kx ■ cos ky pairing state are qual- 
itatively similar to those in the usual s-wave state. How- 
ever, the sign-changed pairing state can sustain more ro- 
bust bound state solutions without a tt phase shift of the 
SC gap near the impurity in the strong scattering regime. 
The rich spectral features in our calculated energy- and 
space-dependent local density of states (LDOS) may be 
directly resolvable by future scanning tunnel microscope 
and assist to ultimately determine the phase structure of 
the SC order parameter. 

The organization of this paper is as follows. In Sec. 
II we briefly introduce the model and formalism we 
adopted. In Sees. Ill and IV we present our results for 
the cases of non-magnetic- and magnetic-impurity sys- 
tems, respectively. Finally, some remarks are given in 
Sec. V and we conclude in Sec. VI. 



II. MODEL AND FORMALISM 

Our tight-binding microscopic Hamiltonian for the 
iron-based superconductors, describing iron atoms on 
a two-dimensional square lattice with two orbitals per 
site, is based on the so-called two-orbital exchange cou- 
pling model developed in Refs. [tI IssIIs^ . Explicitly, 



Hq = Hi + Hint , where the non- interacting part reads 

Ht = ^vl/t(k)/,^(k)vl;(k), 
k 

ht{k) = [(e+(k) - ^)o-o + e-(k)cr3 e^j^(k)CTi] (g) T3 

(1) 

with 'l'^(k) = (c|_^ ^, ci,_k,i, 4 |, C2,_k,i) in Nambu 
spinor representation, j, ^ creates an electron car- 
rying momentum k with orbital a {a = 1,2 for two 
degenerate "dxz^' and "dyz" orbitals, respectively) and 
spin polarization a. For a compact notation, we have 
made use of two sets of Pauli matrices, cr^ and t^, act- 
ing on orbital and particle- hole spaces, respectively, with 
(To or To the 2x2 identity matrix. The matrix el- 
ements of ht, e+(k) = — (ti + t2){coskj; + cos ky) — 
cos kx cos ky, e_(k) = — (t^ — t2)(cos kx — cos ky), and 
£xyO^) = — 4t4 sin kx sin ky are parametrized by four hop- 
ping amplitudes. The normal-state Fermi surfaces in 
the unfolded BZ can be reasonably produced by setting 
ti = -1.0, t2 = 1-3, andta = <4 = -0.85 (see Fig. [T]). For 
convenience, we have taken = 1 as our energy units, 
lattice constant a = 1, and also have made fj, = 1.65, 
which corresponds to electron density ~ 2.1 per site. 
The interacting part contains several terms and can be 
expressed as 

Hint — ^ ^ ^ ^ '^l(SQ^r ' Sa,r' ^^a,r^a,r') 
(rr') a 

+ ^ ^ ^ ] ■^2(Sct,r ■ Sc^r' — ^'^a,r^a,r') 
((rr')) " 

+ (2) 

where Sa.r ~ af^aa'Ca.r.a' and n^.r arc the local 

spin and density operators with orbital a. (rr') and 
{{yy')) denote nearest-neighbor (NN) and NNN pairs of 
sites, respectively. "• • • " represent our ignored orbital 
crossing exchange coupling and Hund's coupling terms, 
which are argued by one of us in Ref. to be unim- 
portant on determining the pairing symmetry of the SC 
state in this model. 

In this paper, we will focus on the s-wave cos kx ■ cos ky 
pairing symmetry in iron-based superconductors and ne- 
glect uncompetitive/subleading pairing symmetries such 
as s-wave cos kx + cos ky and d-wave cos kx — cos ky (as 
shown in Ref. 0) by setting exchange coupling = 0. 
Also, we will assume that the low-energy physics of the 
system may reliably be captured by the mean-field ap- 
proximation as long as the pairing interaction is small 
compared to the bandwidth. By defining the local s±- 
wave pairing amplitude for orbital a, 

AQ(r, Y + 5) = - J2(Ca,r,iCa,r+5,T) (3) 

with 5 = ± y indicating NNN pairing, the mean-field 
Hamiltonian of i^o is then written as 

i?o^^^ = i/t-^i 5^[A;(r,r + 5)c„,rac„,r+5.T + /i.c.]. (4) 
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FIG. 1: Fermi surfaces of the two-orbital model representing 
iron pnictides at = 1.65 in the (unfolded) FBZ. The blue 
(red) curves correspond to the hole (electron) Fermi pockets. 
The ± sign within pockets indicate the relative sign change 
in the cosfc^, -cosfcy SC order parameter. Also, the solid thin 
lines mark the nodal lines of it. 



Furthermore, as a useful comparison, we shall consider 
the "sign-unchanged" onsite s-wave symmetry as well in 
order to extract the non-triviality of the "sign-changed" 
s-wave symmetry. In this case, we replace the interacting 
term by Hint = -\V\Yl,a,v'^a,T,\na^r,i and make use of 
the following mean- field decoupling instead. 



A. Self-consistent Bogoliubov-de Gennes 
formulation 

In order to treat the impurity scattering problem and 
a spatial variation in the superconducting order param- 
eter on equal footing, we will mainly use self-consistent 
BdG formulation,—!^ which has been proved beneficial 
in gaining real space information, to demonstrate our re- 
sults of investigation. 

Within this formulation, we diagonalize the quadratic, 
mean-field Hamiltonian ^ plus impurity term [either 
Eq. (O or dH])] through the BdG equation, 



A| 

V 





A2 



I 



n 



(9) 

with nth eigenvalue and the operators in the matrix 
above obey 



+ [{Vi +aJis,/2)S, 



t2'"2r±xcr il"2r±j)<T 



- ^3 51 <r+Sa 



+ [{Vl + aJ7S./2)(5r,r, - A«]«2r., 



''a.r.cr' ''Ct.r 



where 



ifl2aMlr, 



*2r+5(T 



{V{+<jj'jS,/2)Sr^rrU"2 



i ^ A„(r, r + S)v2r+sa^ (similar to 



Asair) = -|C/|(Ca,raCQ,r,T), 



(6) 



The interaction between the conduction electrons in 
the superconductor and a single non-magnetic impurity 
on site r/ can be written as 



(7) 

where Vi (V/) represents the intra-orbital (inter-orbital) 
scattering strength. On the other hand, the scattering 
from a static (classical) magnetic impurity with magnetic 
moment s is given by 



H. 



mimp — I ^ ^ 



(8) 

where J/ (J|) represents the intra-orbital (inter-orbital) 
magnetic scattering strength. Note that due to spin- 
rotational symmetry, one can choose the coordinate sys- 
tem for the spin degrees of freedom such that z axis points 
in the direction of s. Consequently, it is sufficient to keep 
only z-componcnt in Eq. ^ hereafter. 



where cr = ± correspond to spin up/down, S are NNN 
vectors, and Q — {^,^). The relation between quasi- 
particle operators 7 and electron operators is Ca^r,tT ~ 

J^ni'^arala.n^a -crv2*^-fi^n,s) , and heuce Combining with 
the definition of s-wave cos fc^; - cos ky SC order parameter, 
this gives rise to the following self-consistent conditions. 



Aa(r,r + ,5) 



J: 



2 V^/ 



X tanh 



2 / .(^arT^ar+aj. + "ar.|.^ar+.5T 
n 

En 



2kBT 

For onsite s-wave pairing, we instead have 



(10) 



E„ 



2^ Para I 



n 



'f{En), 
(11) 



where f{E) is the Fermi distribution function. No- 
tice that the summation here is only over those eigen- 
states with positive eigenvalues due to the symme- 
try property of the BdG equation in the whole spin 
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space: 



If 



f)* is an eigen- 



Em then 
s also an 



function of the equation with eigenvalue 

(,,n* 11* n* n* n* n* n* „.n*\t 
1 "21 ' "li I "2^1 '-ll' '-'24 I "it ' "2t / 

eigenfunction with eigenvalue —En- 

Unless otherwise stated, we always perform our com- 
putations on a finite lattice of = 32 x 32 sites with pe- 
riodic boundary conditions at temperature fc^T = 0.03. 
We obtain the resulting quasi-particle spectrum by re- 
peatedly diagonahzing BdG Eq. ^ and iteration of the 
pairing amplitudes according to self-consistency condi- 
tion PH]) [or Eq. pT|) for onsite s-wave] until sufficient 
accuracy is achieved {e.g. the relative error of the pairing 
amplitudes is less than 1%). We choose suitable J2 = 8 
and U = 2.56 such that the ratio of the SC coherence 
peak, Acoh ~ 0.4, to the maximum bandwidth of the 
bands, Wmax = 12, is around 0.033. With this choice, we 
can restrict the SC coherence length ^ ^ hvp / Acoh S 4a, 
as suggested in experiments. Finally, it is practically use- 
ful to note that for the case of non-magnetic impurity, 
the subindex a becomes unimportant for u and v, and 
we can save the computation effort by cutting the spin 
space half, i.e., keeping only index a =t in Eq. 
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B. Non-self-consistent T- matrix approximation 

By assuming that the relaxation of the superconduct- 
ing order parameter is negligible, we greatly simplify our 
impurity-scattering problem and make the physics more 
transparent. As we will see later, the validity of us- 
ing T-matrix approach^^*^ to capture qualitatively right 
physics is justified by comparing with the results from 
the BdG formulation. In addition, the main difference 
between a "self-consistent" treatment and a "non-self- 
consistent" treatment in this approach is the inclusion of 
(normal/anomalous) self-energy corrections to the bare 
electron Greens function or not. Basically, these correc- 
tions are proportional to the impurity concentration and 
hence can be reasonably ignored for a single-impurity 
problem, as in our case here. However, the price we have 
to pay for this simplification is the loss of complete in- 
formation around the impurity site. 

Let us start with a mean-field Hamiltonian in the 
SC state, H^^^ = *^(k)/i(k)«'(k), where h{k) = 
/it(k) -I- A(k)cro ® Ti, with notations defined before in the 
momentum space. Note that we have taken a suitable 
gauge choice to make SC order parameter real and set 
for s±-wave pairing, AQ(k) = A(k) = Aq cos fc^, cos fcy, 
while for onsite s-wave pairing, AsQ(k) ~ A(k) = Ago. 
In the same compact notation, now the impurity poten- 
tial can be unified as i/,;,„p = X^k k' ^'^(k)F^(k') with 
V = V^vGfj, ®Ti, (no summation over ^ and v), where dif- 
ferent types of impurity scattering problems are related 
by V^o3 = Vi,Vi:, = l^, I/qo - J7S,/2, and - J|sj2 
(otherwise, V^i, = 0). 

Defining Cj = uj + iO+, the Green's function for a clean 
SC system is given by 



G"(k,t;)) = [Cbh - h{^)\ 





(0/lt,l 



(12) 

where is a 2 x 2 matrix acting on the particle-hole 
space. The full Green's function in the single-impurity 
problem within T-matrix approximation is then written 
as 



G(k,k',t:;) = G"(k,tD)(5k,k'+G°(k,t:;)r(k,k',c:;)G°(k',(i) 

with the whole impurity-induced effect contained only in 
the T-matrix. Standard perturbation theory gives 



FIG. 2: (Color online) LDOS (red curve) at the nonmagnetic 
impurity site, vi = (16, 16), for (a) the s±-wave pairing state 
and (b) the onsite s-wave pairing state with Vi = 4. The 
dotted black curves represent the LDOS at the same position 
(renormalized by a factor of 4) in a clean system for compar- 
ison. (Acoh ~ 0.4) 



T{Cb) = V + Vg^{Cj)V + Vg^{Cb)Vg^{Cj)V + . . . 

= [h-Vg'^m-^V, (13) 

where g"{u}) = / ^|^G°(k, (D). These equations allow 
us to determine the solutions of the impurity-induced 
bound states via Dct(T~^)=0 in the sub-gap regime, 
H < |A(k)|. 
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III. NONMAGNETIC IMPURITY 

We begin with our discussion on the effect of nonmag- 
netic (scalar) impurity in s±-wave superconductors. It is 
weh known that nonmagnetic impurities in a single-band, 
isotropic s-wave superconductor are not hard pair break- 
ers and hence the impurity-induced spectral feature lies 
essentially at the gap edge.— i^i^ However, as we will see 
later, it is not the case for a s±-wave superconductor due 
to its non-trivial SC phase structure in the momentum 
space. Any potential scattering between hole and elec- 
tron Fermi pockets may destroy the phase assignment 
and leads to the formation of nontrivial in-gap bound 
states. 




(O/ltl 



(b) diagonal direction 




co/lt,l 



FIG. 3: (Color online) LDOS as a function of energy uj and 
the distance R away from the nonmagnetic impurity along 
(a) -\-y direction and (b) diagonal direction. The red, pink, 
and green curves represent LDOS at the impurity site, its 
first neighbor, and its second neighbor along y or diagonal 
direction, respectively. All curves with different _Rs are shifted 
by 0.05 along vertical axis with each other. Note that the 
parameters are the same as those used in Fig. [2] and the self- 
consistent pairing potential at r/ is around 0.17. 



Consider first the case of intra-orbital potential scatter- 
ing, where V/ ^ and V/ = J/ = J| = 0. For a s±-wave 
superconductor, we show in Fig. [Slja) the LDOS, which 
is obtained via 

NA^.V) = -l^[|<rjV'(i^n-C.) + |<r.P/'(^^n+C.)], 

(14) 

then summing over spin tr at the impurity site, and com- 
pare it with the one without any impurity (dotted curve). 
For the demonstration purpose, we choose a moderate 
scattering strength, V/ = 4. The spectroscopic signa- 
ture of bound-state solutions is clearly seen as two peaks, 
symmetric with respect to zero energy and within the SC 
coherence peak Acoh ~ 0.4, in the LDOS at the impurity 
site. Furthermore, the weaker spectral weight at positive 
energy and the stronger one at negative energy arise from 
the absence of particle- hole symmetry in the system. The 
presence of such in-gap bound states in s± pairing state 
is indeed in sharp contrast to the similar problem in the 
sign-unchanged s-wave pairing state, where no peak can 
be found within the SC gap as shown in Fig. ^b) . The 
localized nature of the impurity-induced states within the 
gap is further proved by showing the LDOS as a function 
of energy u) and the distance i? = — r/| off the impurity 
position along certain directions in Fig. O As i? is away 
from the impurity site by one or two lattice constants, 
the peaks disappear quickly and the LDOS recovers back 
to the shape of the bulk DOS. 




FIG. 4: Bound-state energy as a function of non-magnetic, 
intra-orbital impurity scattering strength: (a) V/ > and (b) 
Vi < 0. Note that the x axis is drawn in logarithmic scale. 
(A,„^ = 0.4) 



Above results can be qualitatively understood by (non- 
self-consistent) T-matrix approach. Combining (7" (a) ) for 
the s±-wave pairing state derived in Appendix and the 
inverse of the scattering matrix V~^, the inverse of the 
T-matrix in Eq. p3p is given by 



=To-7oT3)], (15) 



where we have defined the dimensionless scattering 
strength c„ = Viirpo [irpQ ~ 0{1) in the iron pnictides, 
and 7o < is related to the used energy cutoff here (see 
Appendix)]. It is important to realize that as < A^^^, 
Im(T^^) 0. Thus, true bound states at real uj could 
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be found by the condition Dct(T ^) 
bound-state energy, 



fin = ±A 



coh 



Cn70| 



•i7oj 



0, leading to the 



(16) 



This is in contrast to the case of the nodal c?-wavc pairing, 
where we usually get virtualhound states at complex 
Note that for each solution in Eq. , it is doubly degen- 
erate due to orbital degeneracy. Also, as sharply opposed 
to the s±-wave pairing, there are no in-gap bound states 
found for the onsite s-wave state, quite consistent with 
the results shown in Figs, ^a) and Hfb). In fact, this 
distinct feature comes from the (nearly) absence of n 
component in T~^(a)), which is proportional to Acoh in 
the sign-unchanged s-wave pairing. 





FIG. 5: (Color online) The spatial dependence of the SC gaps 
for (a) the s±-wave pairing symmetry and (b) the onsite s- 
wave pairing symmetry in the presence of a non-magnetic 
impurity at r/ = (16,16). 

In addition, by increasing the scattering strength 
Vi > 0, the bound-state solutions for the s±-wavc pair- 
ing state, as estimated by T-matrix approach, evolves 
from the g ap edg e in the weak scattering limit to 
±Aco?i7o/\/7o + 1 in the unitary scattering limit. How- 
ever, the qualitative behavior changes greatly when 
increasing the magnitude of the negative scattering 
strength. From Eq. pB]) . reversing the sign of c„ tells 
us that there must be a minimum bound state energy, 
occurring at a critical Vi < 0. (It does not reach zero 



energy in our system because a small ti component in 
would appear practically due to imperfect cancella- 
tion of the A terms on electron and hole Fermi pockets 
and the SC gap relaxation should be taken into account 
as well.) The numerical results, as seen in Figs. IH^a) 
and m^b) , indeed follow what we have discussed from T- 
matrix consideration. 




-1.0 -0.5 




FIG. 6: (Color online) LDOS as a function of energy uj (red 
curves) in the presence of a magnetic impurity at rj = (16, 16) 
for (a) s±-wave with J/Sz/2 = 2 and (b) onsite s-wave with 
JiSz/2 — 2.2. The black, dotted curves represent the bulk 
DOS normalized by the number of sites for useful compar- 
ison (Acoh ~ 0.4). 

Next, we consider the change in the SC gap function 
caused by a single nonmagnetic-impurity scattering. Fig- 
ure Efa) shows the self-consistent SC pairing potential 
A(i") = s ^a{j^,J^ + S)/8 with s±-wave symmetry. 
Since the impurity potential is short-ranged, the SC gap 
changes largely in the vicinity of the impurity site and 
recovers soon to its maximum value away from the impu- 
rity. It is clear that there are two relevant length scales 
controlling the behavior. The shorter one associates with 
the range of the impurity potential, in which the SC gap 
is strongly suppressed due to much smaller electron popu- 
lation. The other scale is comparable to the SC coherence 
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length ^ < 4a, in which the SC gap is weakly oscillating. 
This oscillation simply indicates the competition between 
the impurity potential and the SC pairing potential. An- 
other subtle feature in Fig. [5)^ a) is that the contour line 
of A is anisotropic and in roughly diamond shape round 
the impurity site. This should be due to the fact that the 
DOS from d^z and dyz orbitals is most likely dominated 
by the (elliptic) electron pockets around (0, tt) and (tt, 0) 
in this model, causing the anisotropy of the SC coher- 
ence length. This feature may only slightly depend on 
the specific form of the pairing symmetry since similar 
anisotropy of the spatial dependence of the onsite s-wave 
SC gap function, As(r) = J^a ^sa(r)/2, is also seen in 
Fig. El^b) . Interestingly, we also observe that the spatial 
distribution of LDOS at fio, N{no,r), respect the same 
anisotropy (not shown). 

As a final remark in this section, we comment on 
the case when there exists small component of the 
inter-orbital, nonmagnetic-impurity scattering, i.e., V = 
(V/ (To -I- V}'cti) (g) T3. We sketch the analysis briefly below 
by using the T-matrix approximation. The easiest way to 
consider this problem is to make a unitary transformation 
in the orbital space such that the transformed scattering 
matrix becomes {Viao + Vla3)(^T3. Also, the transformed 
5° ((D) is the same as before due to the fact that it is diag- 
onal in the orbital space. Thus, we now simply deal with 
new impurity potentials Vj ± separately within each 
(transformed) orbital. The direct consequence is simply 
the breakdown of the orbital degeneracy such that each 
bound-state energy f2o splits into two. This is still a 
distinguishable feature from the s-wave pairing state, in 
which essentially no sharp peaks in LDOS at the impurity 
site when \uj\ < Acoh- 



IV. (CLASSICAL) MAGNETIC IMPURITY 




a)=-0.05|tj 




(B=-0.05|tJ 



FIG. 7: (Color online) The spatial distribution of LDOS for 
(a) s±-wave at Qq = —0.05, and (b) onsite s-wave at Qq = 
—0.05. They correspond to the arrows indicated in Figs. EJa) 
and[6l^b), respectively. 



We now turn to the discussion on the effect of magnetic 
impurities in the s± pairing state. The magnetic impuri- 
ties, the pair breakers of the Cooper pairs, are known to 
induce in-gap bound states in the conventional nodclcss 
superconductorsj ^^i^^'^° These bound-state solutions are 
usually localized near the impurity, possibly with non- 
trivial spin configuration around it, and may dramati- 
cally modify the ground-state properties of the system 
once the magnetic interaction is much stronger than the 
condensation energy, Acoh- 

To investigate possible features in our target pairing 
state, we consider the case of intra-orbital, purely mag- 
netic impurity scattering, where Jj ^ while Vj — Vj = 
J'j = 0. Examining the LDOS spectrum at the impurity 
site in Fig. Elja), one can immediately see that there are 
two sharp peaks, symmetric with respect to zero energy 
within the SC gap [Acoh ~ 0.4), indicating the presence 
of the bound-state solutions. The asymmetric spectral 
weights for the peaks are again due to the breakdown of 
the particle-hole symmetry in the system. However, the 
presence of in-gap peaks in the LDOS is also observed in 



the onsite s-wave pairing state as shown in Fig. [UJb) , ex- 
cept for the reversed magnitudes of the spectral weights 
on the two peaks. This subtle difference may not be con- 
sidered as a general feature since the asymmetry of the 
spectral weights depends on the position and the scat- 
tering strength. Furthermore, as we examine the LDOS 
spectra away from the impurity, the peak positions do 
not change and their spectral weights decay rapidly after 
a few lattice constants comparable to the SC coherence 
length ^. In Fig. [7]Ja), the spatial distribution of LDOS 
at fixed bound-state energy, A^(±J7o,r) [see Eq. (fH]) ]. 
further confirms our observation. While in Fig.[7jb), we 
find no essential difference for the sign-unchanged s-wave 
pairing state, except that the spectral weights are more 
concentrated on the impurity site. For energies outside 
the SC gap, the spatial distribution becomes much more 
extended with only a suppressed region around the im- 
purity (not shown here), as one can expect. 

Similar to the nonmagnetic-impurity problem dis- 
cussed in the last section, the T-matrix approach may 



FIG. 8: (Color online) The energy of bound states fio as a 
function of the effective magnetic impurity moment J/Sz/2 for 
(a) s±-wave pairing and (b) onsite s-wave pairing. (Acoj, ~ 
0.4) 



JiSz/2=1.5|tJ 



assist us to understand the properties of the magnetic- 
induced bound states more clearly. For the s±-wave pair- 
ing symmetry, the inverse of the T-matrix is obtained via 
the replacement of Vf^ and c„ in Eq. by the inverse 
of the effective magnetic moment, (s'^f^)'^^, and the di- 
mcnsionlcss magnetic scattering strength, = 
(s^fJ^ = J/s^/2). Det(r-i) = is satisfied at^ 



^coh 



1 T |cm7o| 



\/(lT|c™7o|)2 



(17) 



For Cm > 0, 7o < 0, we will ignore the second bound-state 
solution above sign in Eq. ((T7)) ] since it is very close 
the gap edge as Cm <C 1, while its magnitude approaches 
to that of the first solution as Cm ^ 1- Therefore, it is 
hardly discernible in our numerical LDOS results. The 
first bound-state solution in Eq. p7|) is to be compared 
with the one obtained from the case with onsite s-wave 
symmetry, where the in-gap bound state occurs at 



l-d(l+7o') 



V[l-c^(l+7o')P+4c2, 



(18) 



Now, we see that in both pairing symmetries, the pres- 
ence of the in-gap bound states is consistent with our 
numerical results shown in Figs. [6lja) and[6l^b). In addi- 
tion, one should notice that, first, each solution is dou- 
bly degenerate due to orbital degeneracy; second, al- 
though the quasi-particle bound-state energy f2o appears 
not symmetric with respect to zero energy, the resonance 
peaks in LDOS are indeed symmetric since each quasi- 
particle state has its particle and hole components at 
— |f]o| and 1 separately 4i So far, in some sense, the 
magnetic-impurity effect could not be a good probe to 
distinguish the sign-changed s-wave symmetry from the 
sign-unchanged one. Nevertheless, there are still "subtle" 
features possessed only by the s±-wave pairing state, as 
we will explain next. 

The in-gap bound state is more robust for the s±- 
wave pairing symmetry in the strong impurity scattering 
regime. In single-band s-wave superconductors, one of 
remarkable properties due to a localized magnetic im- 
purity is that the first-order quantum phase transition 




JiSz/2=4|tJ 



FIG. 9: (Color online) The spatial distribution of the spin- 
resolved LDOS for s±-wave pairing symmetry at uj = \Q,o\ 
with (a) JiSz/2 — 1.5 before the transition, and (b) JiSz/2 = 
4 after the transition. Only spin-down (spin-up) distribution 
is shown in (a) [(b)], while the others are nearly zero on the 
entire lattice. 



takes place as the effective moment, s'^^f = JiSz/2, is 
greater than certain critical value, s'^^-^ . This transition 
represents the jump of the spin quantum number of the 
ground state from 0, where the magnetic impurity is un- 
screened, to spin 1/2, where the magnetic impurity is par- 
tially screened. We refer interested readers to the review 
paper by Balatsky et al. for detailed discussions4^ Here, 
in our case, one indication for such a transition is given by 
Eq. (1171) . where iln switches sign as > — \ — r. This 
is similar to the case of s-wave state but with slightly 
different critical value. However, in the strong scattering 
limit, where Cm 3> 1, we observe that the bound-state 
energy for s±-wave pairing never evolves back to the gap 
edge, but it does for the onsite s-wave pairing, as shown 
clearly in Figs. [D^ a) andlljb), respectively, from our nu- 
merical study (flipping the sign of Cm does not change the 
result). This feature may suggest the bound state solu- 
tion in the sign-changed SC state is more robust than the 
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one in the sign-unchanged case in the strong scattering 
regime. Note that the non-crossing to the zero energy 
is understood due to the local SC gap relaxation and 
many-body effect 



(a) 




but not for the s±-wavc pairing state in Fig. [TlTa). 



Jis=/2 





(b) 


■ 
■ 




1 


s 




1 — 1 


JiSi/2 



FIG. 11: The SC order parameter (dxz orbital only) at the 
impurity site r/ as a function of the effective magnetic mo- 
ment JjSz/2 for (a) s±-wave and (b) onsite s-wave. Note that 
X axis is in logarithmic scale and the rescaling factor Ao of 
the SC gap is chosen to be 0.54 (0.38) for graph (a) [(b)]. 

Finally, we comment on the case when there exists 
small component of the inter-orbital impurity scattering, 
i.e., Jj ^ and Jj ^ 0. Similar to the case with a 
nonmagnetic impurity, this would lead to the splitting 
of the bound-state energy due to the breakdown of the 
orbital degeneracy. Consequently, we now have two crit- 
ical values, s'^(^ and s^g'^, when increasing the effective 
magnetic moment, JiSz/2 with fixed ratio Jj/Jj. These 
transitions correspond to the change in the spin quantum 
number in the ground state, that is, from to 1/2 and 
1/2 to 1, respectively. When sl(-^ < s*^^^ < s]JJ ^ the spin 
and orbital degrees of freedom are strongly correlated and 
more complicated spin configurations are expected. This 
may deserve further study in the future. 





FIG. 10: (Color online) The spatial dependence of the SC 
gaps for (a) s±-wave pairing symmetry and (b) the onsite s- 
wave pairing symmetry. The parameters are same as those 
used in Fig. [S] 

A further evidence to verify this transition (level cross- 
ing) can be provided by the spin-rcsolvcd LDOS, de- 
fined in Eq. (|14p . around the impurity site. As shown in 
Figs.IHa) and[9jb), one can clearly see that the spin con- 
figurations are completely interchanged as JiSz/'i crosses 
the critical value~ 2.2. Similar results are also found for 
the onsite s-wave pairing state (not shown). 

No IT phase shift of the SC gap function around the im- 
purity for s^^-^ > s'^J^ with the s±-wave pairing symme- 
try. In Figs. [TOT a) and fTUT b). we show the self-consistent 
SC pairing potentials for the s±-wavc and s-wave pair- 
ing symmetries, respectively. The basic features arc not 
much different from the cases we have discussed in the 
non-magnetic impurity problem. However, as pointed 
out by Salkola et al.^ the phase of the SC order pa- 
rameter changes by tt at the magnetic-'mrpwcity site with 
respect to the bulk phase when s'^^^ is larger than the 
critical value. This is indeed a sharp feature we have seen 
for the sign-unchanged s-wave pairing state in Fig. lllf b). 



V. REMARKS ON NODAL PAIRING 
SYMMETRIES AND MANY IMPURITIES 

Although wc have mainly concentrated on the effects of 
a single nonmagnetic (magnetic) impurity in fully gapped^ 
sign-changed, s-wave superconductors, we would like to 
make two remarks which are closely related to, or sug- 
gested by the current study. 

(1) Considering the case of sign changed, but with SC 
nodes on electron Fermi pockets, as called for attention 
in recent theoretical works^^i^i^ is of course beyond the 
scope of this paper and may deserve a future study. How- 
ever, following the same T-matrix calculations sketched 
above, wc can consider certain nodal pairing states in a 
straightforward manner. For instance, the gap function 
could be of the form. AQ(k) — Ao(cosA:j; -f cosfcj,)/2, 
i.e., the Aig nodal pairing symmetry with nodes on 
the electron pockets. When evaluating (?°(u)), one can 
realize that the ri component in Eq. (jl5p is not ab- 
sent (or negligible) anymore (with contribution mainly 
from the hole pockets), and it makes the system much 
closer to a sign-unchanged s-wave pairing state. Thus, 
in the single nonmagnetic-impurity problem, the should- 



10 



bc-prcscnt peaks in the LDOS disappear or are nearly 
indistinguishable with the gap edge (continuum). 

(2) There are already several paper o^^i^^'^^ discussing 
about the issue of many impurities and its relation to 
the SC Tc with sign-changed s-wave pairing. In partic- 
ular, Senga and Kontani^ present a detailed study on 
the (nonmagnetic) impurity-induced DOS and the sup- 
pressed Tc with various inter-band (intra-band) scattering 
strengths /' (/) within T-matrix approach. According to 
their results, at the fixed ratio I' /I = 1, the induced- 
impurity band would move toward zero energy without 
going back to the gap edge as increasing the scattering 
strength, associated as well with a large suppression of Tc- 
This tendency corresponds to our observation that the 
impurity-induced bound-state energy never evolves back 
to the gap edge as increasing the scattering strength, and 
should be in sharp contrast to the sign-unchanged s-wave 
pairing state, in which Anderson's theorem^ is expected 
to be satisfied. 



VI. CONCLUSIONS 

In conclusion, we have studied the impurity-induced 
in-gap bound states in iron-based superconductors with 
(sign-changed) s±-wave pairing symmetry by using 
both the self-consistent BdG formulation and non-sclf- 
consistent T-matrix approach. In comparison to the sign- 
unchanged s-wave pairing state, we have found several 
signatures, which are mainly associated with the sign 
change in the SC order parameter. In particular, for a 
nonmagnetic impurity, the two in-gap bound-state peaks 
appear in the LDOS at or near the impurity site and 
their formation is due to the sign-reversal effect in the or- 
der parameter during Andrecv reflection processes. For 
a magnetic impurity, there also exist bound-state solu- 
tions, but only for one of the electron-spin polarizations 
around the impurity at the resonance energy due to the 
breakdown of the local time-reversal symmetry. Above 
a critical effective magnetic moment, the ground state of 
the system undergoes a quantum phase transition, from a 
spin-unpolarized state to a spin-polarized one. Although 
in the presence of a magnetic impurity, both sign-changed 
and sign-unchanged s-wave pairing states behave quali- 
tatively the same, we emphasize that the former pairing 
state can sustain more robust bound-state solutions with- 
out a TT phase shift of the SC gap near the impurity in 
the strong scattering regime. 

Note added- As we nearly complete our paper, we find 
two interesting papers by ZhangylS, and Tao Zhou et al.fi^ 
discussing similar issues on the impurity effect. Their 
insightful results are basically consistent with ours. 
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APPENDIX A: DERIVATION OF g°{Cj) WITHIN 
T-MATRIX APPROXIMATION 

In order to calculate the bare Green's function G°(k, uj) 
and its corresponding g^[Cj), it is more convenient to turn 
our orbital basis into band representation, where we can 
easily obtain the Green's function for each band, by the 
following unitary transformation 



^t(k)( e (k) 



where 



C/(k) 



cos dk 



e^y(k) e^y(k) 



cos(6ik/2) -sin(6lk/2) 
sin(0k/2) cos(0k/2) 

e_(k) 



^e(k) \ 

e/.(k) ) 



(A2) 



6^(k)+62 (k)' 



= , sm ffe 



e^y(k) 



62„(k)+e2 (k) 



and e^(j,)(k) = e+(k) ± e_(k). In the normal state, 
ee(ft)(k/) — /i associates with two electron (hole) pockets. 
H^'P now transforms as H^^ = Xik *^(k)/j(k)*(k), 
where /i(k) = {[(ee(k) + e,,(k))/2 - ^Ji]ao + [(ee(k) - 
e/i(k))/2]cr3} (E) T3 -f A(k)o-o ® ti. Consequently, the bare 
Green's function G'^(k, a)) for band electrons is given by 

[(. + .0-)/.-i(kr..(«S(J-)g.;^_^,),,A3) 

where 

n ~\ '^7-o + A(k)Ti + (ec(^)(k)-^)T3 

Ge(h) k, w = — — — \^ . A4 

^ ' ^~ (ee(/o(k) - W - A^(k) 

Now, we can transform above Green's function back to 
its orbital representation defined in Eq. (|12[) with matrix 
elements, 



G?i(22)(k,'^) - G°(,)(k,^)cos2(Ji) 
-I- G$l(,)(k,i;)sin2(^), 



(A5) 

G;2(2i)(k,<i) = [G°(k,(i)-G",(k,cD)]sin(^)cos(^). 

Upon integrating over momentum in the FBZ, we ob- 
tain the following approximate expression for g'^ioj), 



d^k f GUK^) GUK^ 



(27r)2 \GUK^) G^Klu 

1 r d^k 

2./ (2^ 



or -npoao > 



a{uj) 



- 7((i)T3 , (s) 



a{uj){ujTo -f AqTi) 



(A6) 
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FIG. 12: The arrows show the direction and magnitude of 
the vector {exy(k),e-(k)) in the (unfolded) FBZ. It is then 
easy to extract 9u, information from this arrow map. The 
orbital-resolved Fermi pockets are also put on top of the map 
for convenience. Red and green colors represent d^z and dyz 
orbitals, respectively, {fi — 1.65) 



where a{Q) ~ ^ 
1,7('^)= 2^1n 



[tan-i 



tan 



-1 ( E, 

with En 



g|+A^-!i!'i ) '^'0' "^ith Kg representing 
the average energy difference between the pocket center 
and the Fermi level and Ec representing the energy cutoff 
with respect to the Fermi level.— po is the density of 
states at the Fermi level. 



Note that to get the second equality, we have taken 
into account the features of 6'k in the FBZ (sec Fig. [T^]); 
to get the third/fourth equality, we have made several 
approximations: (i) for each pocket, the energy disper- 
sion is quadratic with respect to the pocket center, (ii) 
the density of states around the Fermi level (and within 
the cutoff energy Ec) for each pocket is a constant poi 
and (iii) for s-t-wave pairing, A(k) ^ Acoh for hole pock- 
ets while — Aco/i for electron pockets. This is in contrast 
to the on-site s-wave case where we take A(k) = Acoh 
for all Fermi pocketSf^^ 
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